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Abstract 

Novel contributions to the total inclusive cross section for Higgs-Strahlung in the 
Standard Model at hadron colliders are evaluated. Although formally of order a^, 
they have not been taken into account in previous NNLO predictions. The terms under 
consideration are induced by Higgs radiation off top-quark loops and thus proportional 
to the top-quark Yukawa coupling. At the Tevatron, their effects to WH production 
are below 1% in the relevant Higgs mass range, while for ZH production, we find 
corrections between about 1% and 2%. At the LHC, the contribution of the newly 
evaluated terms to the cross section is typically of the order of l%-3%. Based on these 
results, we provide updated predictions for the total inclusive Higgs-Strahlung cross 
section at the Tevatron and the LHC. 

1 Introduction 

With the LHC experiments becoming sensitive to signals for the Standard Model (SM) 
Higgs boson, the search for this elusive particle has entered a new and hopefully its final 
phase. The direct searches at LEP, Tevatron, and LHC only leave relatively small allowed 
windows for the Higgs mass Mh, the widest one between 114 GeV and about 145 GeV (see 
Ref. (T|[2] for preliminary results). 

As opposed to evidence or discovery, the exclusion limits rely heavily on theoretical predic- 
tions. The dominant cross section to compare the experimental measurements to is gluon 
fusion which receives large radiative corrections. Although it is probably the most-studied 
cross section for an unconfirmed particle, the residual theoretical uncertainty is still siz- 
able and highly disputed (for a recent discussion, see Ref. |3]). A lot of this uncertainty is 
induced by quantum chromodynamics (QCD), specifically the strong coupling Og and the 
parton density functions (PDFs). 

For various reasons, however, gluon fusion need not be the dominant search mode. The 
focus of this paper is on the associated production of a Higgs boson with an electro-weak 




(a) (b) (c) 

Figure 1: (a) Leading order Feynman diagram contributing to the Higgs-Strahlung pro- 
cess; (b) real corrections at NLO QCD; (c) gg component not covered by Drell-Yan-like 
corrections. 



gauge boson [pp — VH, V G {VF^, Z}), or "Higgs-Strahlung" for short. At the Tevatron, 
where the 77 and 65-decays of a Higgs boson produced in gluon fusion cannot be separated 
from the background with sufficient precision, this mode is particularly important in the 
low mass region. Nominally known through order a^, we identify and evaluate a previously 
neglected contribution which formally adds to these next-to-next-to-leading order (NNLO) 
terms. 

For the LHC, the relevance of the Higgs-Strahlung process used to be considered marginal. 
This has changed with the idea of focusing on events with highly boosted Higgs bosons 
by analyzing the substructure of jets [4]. Even though for a proper theoretical prediction 
in such an analysis one needs to consider differential quantities, it is important to ensure 
that all effects that contribute to the total rate are under control. 

Once the Higgs mass is known, precise predictions for the individual production and 
decay channels will be essential in order to extract the maximum information from the 
experiments (see, e.g., Ref. [5]). 

The theoretical prediction of the total inclusive cross section due to Higgs-Strahlung at 
hadron colliders is under very good controQ the leading order contribution is of order g^, 
where g is the weak coupling constant, and is completely analogous to what used to be 
the main search channel at LEP, except that the initial e~'~e~ is replaced by a qq' pair, of 
course, see Fig.[l](a). The NLO QCD and the bulk of the NNLO QCD corrections |], 
i.e. 0{g^as) and 0{g'^ag), can be reduced to the Drell-Yan production of a virtual gauge 
boson (9|[To]. The theoretical uncertainty due to PDFs has been estimated to be at the 
percent level, and the renormalization/factorization scale dependence of these terms is 
practically negligible For ZH production at O(a^), however, there are a few classes 
of diagrams that have no correspondence to the Drell-Yan process. For example, the 
gluon-induced virtual corrections mediated by a top-quark loop are of order g'^X'^a'^ (see 
Fig.[l](c)), where Aj is the top-quark Yukawa coupling which, in the SM, is of order one. 
These corrections were evaluated in Refs. (8, 12 and found to be of the order of 5% at the 



^For recent work on higher order differential WH cross sections, see Ref. jo]. 
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LHC. 

In this paper we consider another class of diagrams which are formally of order g'^Xta'l 
and were neglected in previous analyses. For simplicity, we will refer to them as "top- 
mediated terms" in this paper, even though they are not the only contributions involving 
top-quarks, as noted above. Their numerical impact is at the percent level and therefore 
within the current estimated theoretical uncertainty of the NNLO result (see Ref. 
Note, however, that this uncertainty estimate is dominated by the effects from PDFs and 
Os ; once these will be known with higher precision, the results of this paper will be required 
for the perturbative part to compete with this precision. 

The remainder of this paper is organized as follows: Section [2] defines the effects to be 
calculated, briefly describes the methods applied, and presents analytical expressions for 
part of the results. In Section [3j we study the size of the newly evaluated effects and 
present updated values for the total inclusive cross section for WH and ZH production at 
the Tevatron and the LHC at collision energies of 7 and 14 TeV. 



2 Calculational details 
2.1 Outline of the problem 

The Feynman diagrams of the top-mediated terms considered in this paper can be divided 
into four groups which will be described in this section. 

Examples of diagrams of the first group, named Vi in what follows, are shown in Fig.[2](a) 
and (b) . They are characterized by the emission of a Higgs boson off a top-quark bubble- 
insertion into an internal (i.e. virtual) gluon line. They contribute to the total cross section 
through the interference with the leading order amplitude (see Fig.[T](a)). 

The second group (-Ri), see Fig.[2](c), can be viewed as the real emission counterpart of 
group Vj. It is obtained by radiating the Higgs off a top-quark bubble-insertion into an 



3 



(a) 



000000 









\ 



.' H 




(b) 




(c) 



Figure 3: (a),(b) Diagrams of group Vu and (c) group contributing to the process 
qq — )• ZH{g) at order g^Xtol. 



external gluon line. These diagrams have to be interfered with the real-emission amplitude 
contributing to the NLO QCD cross section, see Fig.[l](b). Needless to say that the crossed 
amplitudes, where the gluon is in- and a quark or anti-quark is out-going, have to be taken 
into account as well. 

The third group of Feynman diagrams (Vii) is closely related to the gluon-induced contri- 
bution of Fig.[T](c). While the latter enters the total cross section with its square, the new 
contributions treated here are genuine two-loop terms which are to be interfered with the 
tree- level amplitude of Fig.[l](a). Examples are shown in Fig.[3](a) and (b). Finally, the 
fourth group (.Rii) can again be seen as the real-emission counterpart of group Vi\. An 
example is shown in Fig.[3](c). Also here, of course, one can cross the gluon and the quark 
or anti-quark from the final to the initial state and vice versa. 

The amplitude for each of these groups is separately gauge invariant and UV- as well 
as IR-finite, despite the fact that two of them can each be viewed as real and virtual 
correspondences of each other. 

The diagrams of Ri and Rn can be calculated exactly, taking into account the full depen- 
dence on the top, Higgs, and vector boson mass. The amplitude arising from groups Vi 
and Vii, on the other hand, is very challenging to compute and may be even beyond cur- 
rent technology. We therefore follow the established and successful method of asymptotic 
expansions in order to approximate the result in the limit of infinite top-quark mass Mf. 
The validity of this approach will be discussed below. 

In a fully consistent treatment of the contributions described in this section, one needs to 
take into account yet another set of Feynman diagrams. They can be obtained from the 
previous ones by radiating the Higgs boson off the external vector boson instead of the 
top-quark loop, see Fig.|4| for examplej^ If all couplings are replaced by their SM values, 
the resulting amplitudes are of the same perturbative order as the ones discussed in this 
paper; however, they receive an additional suppression factor ~ My/{s — My), where 
My is the mass of the vector boson and \/l > My + Mh the partonic center-of-mass 



^Due to the requirement of anomaly cancellation, however, they have to be supplemented by triangle 
diagrams with bottom- instead of top-quarks running in the loop. 
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energy. Justifiably so, these terms have been considered Drell-Yan-like in Ref. jSj; in fact, 
they can be calculated by convolving the output of (a suitably modified version of) the 
program MASSIVE (l3| with the decay V* — )• VH. As a check, we also calculated some 
of these contributions directly using the methods described below. The result confirms 
the statement of Refs. [8, 12 that these contributions are numerically irrelevant: they are 
typically 2-3 orders of magnitude smaller than the other newly evaluated contributions 
considered in this paper, and will therefore be neglected in what follows. 



2.2 Calculation of Vi and Vu 



A priori, an expansion of the two-loop amplitudes in the limit of large top-quark mass 
seems unjustified, because the partonic center of mass energy y/§ at the Tevatron and the 
LHC can be much larger than Mf (or rather 2Mf which corresponds to the threshold and 
is thus the relevant scale). Nevertheless, there is a number of arguments that make such 
an approach reasonable, if Mh is not too large: 



The parton luminosities at large s are strongly suppressed, and the bulk of the 
contribution to the total cross section indeed arises from the region below the top- 
quark threshold. This is similar to what happens in the case of gluon fusion, see 
Refs. [l4-16 . For Higgs-Strahlung, however, the situation is somewhat worse, be- 



cause the energy window for which the heavy-top expansion is expected to converge 
is much narrower than for gluon fusion: My + Mh < < 2Af(. In fact, for 
Mh ^ 250 GeV, this restriction cannot be obeyed at all. 

For group Vi (Fig.[2](a),(b)), the typical energy scale affecting the top-quark loop is 
not the full center-of-mass energy, but significantly below that, because the vector 
boson carries off a large fraction of the momentum. 

Since group Vu (Fig.[3](a),(b)) is closely related to the 5f(7-induced contribution of 
Fig.[T](c), we may estimate its impact by the ratio of the qq and the gg luminosity 
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Figure 5: Asymptotic expansion of the diagram in Fig.[2](a). The diagrams left of ^ 
are evaluated after setting their external momenta to zero. The result determines the 
expression to be inserted into the effective qqVH or ggH vertex in the diagram right of 

for example. 



For details on the general method, see Refs. 29-31 



times the 5fg'-induced cross section. It should therefore not exceed a few percent of 
the LO cross section. 

Overall, we expect the top-mediated terms to affect the cross section at the percent 
level [6 , 17 . The leading term of an expansion in terms of 1/Mt should thus be 



sufficiently precise. 

The tools to calculate the diagrams are by now standard: asymptotic expansion of the 
two-loop diagrams of groups Vi and Vu leads to a factorization of scales into either two- 
loop massive tadpoles (i.e., vanishing external momenta) times tree-level diagrams with an 
effective qqH or qqVH vertex, or one-loop massive tadpoles times massless one-loop dia- 
grams with an effective ggZH or ggH vertex. As an example, the graphical representation 
of the asymptotic expansion of the diagram in Fig.[2](a) is shown in Fig.jsj 

We use the automatic setup consisting of qgraf [l8| for the generation, and q2e/exp (lOpO 



for the expansion of the diagrams, as well as MATAD ^21j for the calculation of the tadpole 
integrals. In order to evaluate the massless one-loop box and triangle integrals corre- 
sponding to the right-most diagrams in Fig. [5], we supplement this setup by an additional 
routine (based on FORM |22] ) which implements standard Passarino-Veltman reduction f23] 
in algebraic form. The scalar one-loop Feynman integrals are evaluated using the results 
of Ref. [241. 



It turns out that for Vj, the sum of terms involving an effective qqVH or qqH vertex does 
not contribute at leading order in 1/Mt, i-e., the corrections can be calculated by simply 
using an effective ggH vertex and evaluating one-loop diagrams. This observation allowed 
us to obtain V[ in a second, independent calculation by using a generalized version of the 



FORM program FDiag ^25j and the Fortran package FF 26-28j, supplemented by a routine 
to facilitate the tensor reduction of rank-4 tensor 4-point functions and the subsequent 
numerical evaluation of the corresponding tensor coefficients. 
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Contracting with the LO amphtude and summing/averaging over final/initial color and 
spin degrees of freedom, we find for the amplitude of Vj, cf. Fig.[2](a),(b): 




+ {s-Mjj + M^) 



In" ( 4 1 + 2Li 



u 



-2Li2 (1-^1+ 2Li2 ( 1 




where s = [pi + P2)^, t = (pi — PvY ^-i^d u = {pi — ph)'^ are the usual (partonic) 
Mandelstam variables, with the incoming momenta pi,P2, and ph,Pv the momenta of the 
Higgs and the vector boson, respectively. My is the mass of the emitted vector boson, 
and Li2 denotes the di-logarithm. The electro-weak couplings are given by 

Gqq' = i'^q + '^q)^qq' ) = - \Vqq'\'^ , 

^1 1 r -4 1 . 2^ ^1 , f U,C ^ 

^'' = ^2 + 3 1+2/^^^^^' "'^ = ^2 

with the weak mixing angle sin^ 9w = 1 — M^^/M^ and the CKM matrix elements Vqq' 
(see Ref. |32j for the latest numerical values; we set Vqq' = if both q and q' carry the 
same weak isospin charge I3). We have expressed the top-quark Yukawa coupling by the 
tree-level relation 



^ ^ = mJV2Gf (3) 
V ^ 

in Eq. ([T]), and also in Eq. ^ below. Similar to the gluon fusion process, the factor Mt 
cancels against its inverse from the top-loop integration. 

For group Vji (see Fig.[3](a),(b)), only the axial vector part of the Z coupling to fermions 
contributes. We implement its Dirac structure with the help of Levi-Civita symbols 



33 



7m75 ^ ^e^,ap57'^l^l^ , (4) 

and re-write their product in terms of the anti-symmetrized (denoted by square brackets) 
D-dimensional metric tensor, 

£ij.vpa£ - g^^g^ gp g„] ■ [o) 
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Upon asymptotic expansion of the diagrams, the terms corresponding to an effective ggZH 
vertex contribute only at subleading order in 1 /Mj , and the calculation reduces to massive 
2-loop tadpole diagrams. 

The result assumes the simple form 

" - F z t Q / M J 2s - M|f + ^ + 0{ — ) } . (6) 



di 127rs2(s - M|) V vr 




Since the expressions in Eq. (|lj) and ([6j) are free of divergences in the allowed f-region, we 
can numerically integrate them together with the convolutions over the PDFs in order to 
get their contribution to the total inclusive hadronic cross section. The numerical results 
will be presented in the next section. 

2.3 Calculation of _Ri and _Rii 

Being of one-loop order, the diagrams of groups i?i and i?n can be calculated including the 
full dependence on the top-quark, Higgs- and vector-boson mass by means of Passarino- 
Veltman reduction to scalar one-loop functions. In most of the phase-space, the numerical 
evaluation of these function can be performed with the help of the programs FDiag and 



FF (see Sec. 2.2), and independently with FeynArts, FormCalc, and LoopTools [34| (the 
latter of which relies on FF though). The contribution Rn, however, involves momentum 
configurations that are outside FF's capabilities. For those, we use an implementation 



of the one-loop integrals based on Refs. 35 36P Again, since projection with the LO 



amplitude leads to a finite expression, phase-space integration and convolution with PDFs 
can be done fully numerically. 

2.4 Validation of the heavy-top approximation 

This section describes both a consistency check for our calculation of the real-emission 
contributions Ri and Ru, as well as a validity check of our approach to the virtual terms 
Vi and Vii. 



In addition to the exact calculation as described in Section 2.3 we evaluated the real 
emission amplitudes Ri and Ru also by applying asymptotic expansions in the heavy-top 
limit. In these cases, it reduces to a naive Taylor expansion of the integrand before loop 
integration. We can again perform all phase-space integrals numerically. In the gq-channel, 
however, we also integrated over the angular variables of the phase-space analytically in 
D = 4 — 2e space-time dimensions, and found the cancellation of all poles at e = 0. The 
numerical integration over the remaining energy variables is straightforward and leads to 
the same result as the all-numerical method. 

^We thank Stefan Dittmaier for providing us with his private code. 
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Comparing this heavy-top result for Ri and Ru at the partonic level to the exact results 
described above, we indeed observe that heavy-top limit systematically approaches the 
exact result as £/Mt — t- 0, where £ is any of the external energy scales. In particular, we 
observe that, as opposed to the other contributions, the terms Ru vanish as in the 

limit M( ^ oo. This agreement between the exact and the asymptotic approach provides 
a strong and valuable check on our results. 

In order to draw conclusions for the validity of the heavy-top expansion, however, the 
relevant quantity to compare is the hadronic cross section, of course. For that, we find 
that the heavy-top result is roughly within 25%/35% of the full result at the LHC/Tevatron 
in the mass range considered in this paper. 

We expect a similar quality of the heavy-top result for the contributions from Vi and 
Vii. To be conservative, we attribute an additional uncertainty of 30%/50% (relative to 
the central values) to these terms for the LHC/Tevatron prediction of the hadronic cross 
section. 



3 Numerics 

3.1 Size of the top-induced terms 

We present numbers for the Tevatron {pp @ = 1.96 TeV), as well as the LHC {pp @ 
7 TeV and 14 TeV). The hadronic cross section is evaluated by folding the partonic cross 
section with PDFs. Although the corrections evaluated in this paper are not renormalized 
by lower order terms, we think it is appropriate to convolve them with NNLO PDFs. We 
use the central MSTW2008 set, implying as{Mz) = 0.1171. For the physical parameters, 
we assume the values 

Mz = 91.1876 GeV , Mw = 80.398 GeV . (7) 

Since the amplitudes considered in this paper are UV- and IR-finite, they do not depend 
explicitly on the renormalization or factorization scale (//r, /ip). However, these scales 
enter implicitly through the strong coupling and the PDFs. For the numerical analysis, 
we choose 

= {pH + Pvf 

as the central scale for and ^Up, where pu and pv are the 4- momenta of the Higgs and 
the outgoing vector boson. For the Drell-Yan-like terms, equals the invariant mass 
of the intermediate gauge boson. In order to estimate the theoretical uncertainty of our 
results, we vary /xr by a factor of three around this central scale while keeping /xp fixed, 

^V^<m<3v^, IJ'F = \/^, (8) 
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then repeat the analysis after interchanging fip and /Ur, and take the extreme values of 
the cross section as uncertainty band. 

Fig. [6] shows the contribution of the newly evaluated terms to the total inclusive cross 
sections a{pp — )• WH) = a{pp W~^H) + a{pp — )• W~H) and a{pp — )• ZH) at 7TeV and 
14 TeV center-of-mass energy. The contributions from the various groups of diagrams (see 



Section 2.1) are shown separately, each of them with an error band derived from the scale 
variation described above. Also shown is the sum of all contributions. The corresponding 
plots for the Tevatron are shown in Fig. [7} 

The size of the Vi component amounts to about 0.5% of the LO cross sectiorj^ both for 
WH and ZH production, independent of collider type, center-of-mass energy, and Higgs 
boson mass. While at the LHC the Ri terms are typically a little larger than that, and 
increase with Mh and the center-of-mass energy, they always remain below 0.5% at the 
Tevatron. 

Since Vi and Ri are the only non-vanishing contributions for WH production, in this case 
the overall effect of the top-mediated terms evaluated in this paper remains below 1% at 
the Tevatron, and ranges between about 1.1% (for Mh = 100 GeV at 7 TeV) to 2.4% (for 
Mh = 300 GeV at 14 TeV) at the LHC (see Fig.[6](a) and Fig.[7](a)). 

For ZH production, also the groups Vu and Ru have to be taken into account. Being 
suppressed by l/M^ at large Mt, as pointed out above, the terms Ru have only a very 
small numerical impact at the per-mille level. The relative contribution of Vu, on the 
other hand, increases with the Higgs boson mass, but remains below 0.8% at the LHC in 
all of the considered Higgs mass range. At the Tevatron, however, it exceeds Vi and Ri 
above Mh ~ 130 GeV and ranges up to about 1% at Mh = 200 GeV. The overall effect of 
Vi, Ri, Vu, and Ru on ZH production at the Tevatron is therefore between 1% and 2%, 
while it ranges between 1.1% (for Mh = 100 GeV at 7 TeV) and 2.9% (for Mh = 300 GeV 
at 14 TeV) at the LHC (see Fig.[6](b) and Fig.[7](b)). 

In each case, a rough estimate of the uncertainty on the top- mediated terms due to scale 



variation is 20-30%. As discussed in Section 2.4, we also include an estimated uncertainty 



on the Vi and Vu terms of 30%/50% for the LHC/Tevatron, arising from the heavy-top 
limit (this is not included in Figs. [g] and [7]). Considering the fact that the corrections are 
at the percent level, this uncertainty will affect the accuracy of the total inclusive cross 
section by roughly 0.5%. We do not expect the PDF uncertainties on the top-induced 
terms to add significantly to the one of the total cross section and will neglect it in what 
follows. 

*Only central values are considered in this discussion. 
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Figure 6: Contribution of the corrections evaluated in this work to the total 
inclusive cross section a{pp — )• WH + X) (left column) and cr(pp — )• ZH + X) 
(right column) at the LHC with 7TeV (upper row) and 14TeV center-of-mass 
energy (lower row). The effects are shown relative to the leading order cross 
section (Tlq- The various bands show the contributions from Vi, Ri, and, in the 
case of ZH, Vu and i?n- The upper band is the sum of all contributions. The 
width of the band arises from the variation of //r and as described in the main 
text. 
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3.2 Total inclusive cross section 



As shown in the previous section, the effects of the newly evaluated terms are of the order 
of 1-3% of the LO cross section on WH and ZH production. Superficially, this lies within 
previous estimates of the uncertainties of the total inclusive cross section at the LHC [11]. 
However, most of this uncertainty is induced by the PDFs. Once systematical effects 
and theoretical ambiguities will have been settled (for a recent discussion, see Ref. [s]) or 
become irrelevant due to the inclusion of LHC data, one may expect these uncertainties to 
shrink significantly. Then, in order to arrive at the most precise prediction for the process, 
the effects of the top-mediated terms evaluated here will have to be included. 

For this reason, we provide an updated prediction of the total inclusive cross section 
for WH and ZH production at the Tevatron and the LHC at 7TeV and 14TeV in this 
section. The numbers for the LHC are obtained simply by adding the top-mediated terms 



to the predictions of Ref. 11 , which have been obtained using the program vhSnnlo |37|, 



based on the calculation of Refs. f8}|10] and the program zwprod fo], and an electro- weak 
correction factor ^38j. The uncertainties on the numbers given here result from linearly 



adding the previous uncertainties 11 to the ones discussed in Section 3.1 



For the Tevatron, an "official" public agreement on the cross section prediction analogous 
to Ref. 11 is not available. We therefore follow the analysis of Ref. 11 using Tevatron 



parameters, and proceed as above. The electro-weak correction factor we read off from the 



plots and tables of Ref. 38 in this case, despite the fact that the parameters in this paper 



are slightly outdated. The uncertainty induced by this procedure should be negligible, 
however. 
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The final results are presented in Tables [T| [2j and [3] for the LHC at 7 TeV, 14 TeV, and 
the Tevatron, respectively. The top-mediated terms slightly increase all the cross sections, 
but also the theoretical uncertainty. The column labeled "Pert" in the tables contains the 



estimate of the perturbative uncertainty as described in Section 3.1 The column labeled 
"PDF+as" displays the error induced by the PDFs and the input value for Og, as described 
in Ref. fill. 



Finally, Figs.|8] and [9] show the K-factors resulting from our compilation, including NNLO 
QCD, electro-weak effects, and the newly evaluated top-mediated terms of this paper. In 



the LHC case, we also show the previous numbers from Ref. 11 , and for the Tevatron we 
provide the analogous numbers with the same input data. 

These plots clearly show that for the LHC, the top-mediated terms affect the WH pro- 
duction cross section at the order of the estimated perturbative uncertainties. In fact, 
note that for ^/s = 14 TeV and above around Mh = 160 GeV, the uncertainty bands of 
the predictions with and without the top- induced terms hardly overlap. Since the uncer- 
tainty bands of the ZH channel are significantly larger than for WH, mostly due to the 
g^r-induced terms of Fig.[l](c), the relative importance of the top-mediated terms on ZH 
is much smaller than for WH. 

At the Tevatron, the K- factor is typically larger than at the LHC, and the top- mediated 
terms have a rather small impact. Also, the difference between WH and ZH in the K- 
factor is much less pronounced than at the LHC, the reason being again the gfg-induced 
terms which have much smaller impact at the Tevatron. 



4 Conclusions and Outlook 

The effects of an as of yet neglected contribution to the Higgs-Strahlung process at hadron 
colliders have been studied. The relevant Feynman diagrams contain closed top-quark 
loops and are of order g^Xta^. The one-loop real-emission contributions were evaluated 
exactly, while the two-loop virtual terms were evaluated in the heavy-top limit. It was 
argued that this provides a reliable approximation to the exact result. The numerical 
impact of the newly evaluated terms is within the current estimate for the theoretical 
uncertainty, but may well become important once the uncertainty induced by PDFs and 
the strong coupling reduces. 

We plan to include the newly evaluated terms in the publicly available program vhSimlo 
in order to provide a tool for the evaluation of the the full C(a^) prediction for ZH and 
WH production. 

Finally, let us remark that similar contributions exist also for the weak-boson fusion pro- 



cess 139 -42 . Parts of them were evaluated in Ref. 43 
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Mh 

[GcV] 


a{WH) 

[fb] 


Pert 

[%] 


[%] 


a{ZH) 

[ft] 


Pert 

[%] 


PDF±as 
[%] 


100 
105 


1197 
1027 


+1.0 
-0.8 
+0.8 
-1.1 


±3.4 
±3.5 


636.8 
549.8 


+1.5 
-1.5 
+1.7 

-1.9 


±3.4 
±3.7 


110 
115 
120 
125 
130 
135 


883.7 
762.0 
662.7 
578.8 
506.1 
443.7 


+0.7 

-1.0 
+0.8 
-1.1 
+0.8 
-1.1 
+0.6 
-1.2 
+0.8 
-1.1 
+1.1 
-0.8 


±3.8 
±3.9 
±3.4 
±3.5 
±3.5 
±3.4 


476.5 
414.6 
363.3 
319.0 
280.7 
247.9 


+ 1.6 
-1.6 
+ 1.7 
-1.5 
+2.0 
-1.6 
+1.9 
-2.0 
+2.0 
-1.8 
+2.2 
-1.8 


±4.1 
±4.2 
±3.5 
±3.5 
±3.7 
±3.6 


140 
145 


389.9 
344.4 


+1.0 
-0.9 
+0.7 

-1.1 


±3.5 
±3.8 


219.5 
195.2 


+2.0 
-2.0 
+2.3 
-2.2 


±3.7 
±4.0 


150 
155 
160 


303.5 
267.7 
231.9 


+0.9 
-1.1 
+ 1.0 
-1.1 
+ 1.0 
-1.1 


±3.3 
±3.5 
±3.8 


173.2 
154.3 
135.0 


+2.3 
-2.0 
+2.6 
-2.0 
+2.5 
-2.1 


±3.6 
±3.6 
±4.0 


165 


213.2 


+1.0 
-1.1 


±3.6 


124.8 


+2.7 
-2.1 


±4.1 


170 
175 

180 
185 


190.7 
171.0 

154.0 
140.5 


+1.0 

-1.1 

+0.8 
-1.5 
+1.1 
-1.0 
+0.9 
-1.2 


±3.8 
±3.8 
±3.5 
±3.5 


112.0 
100.8 

90.34 
82.46 


+2.7 
-2.3 
+2.7 
-2.3 
+2.8 
-2.3 
+2.9 
-2.4 


±4.2 
±4.1 

±3.8 
±3.8 


190 
195 


126.9 
115.2 


+1.1 
-1.1 
+1.2 

-1.0 


±3.7 
±3.7 


74.65 
67.91 


+2.8 
-2.6 
+2.9 
-2.4 


±3.9 
±4.0 


200 


104.6 


+0.9 
-1.3 


±3.8 


61.81 


+2.9 
-2.4 


±4.1 


210 
220 


86.70 
72.38 


+ 1.0 
-1.1 
+0.9 
-1.3 


±3.7 
±3.7 


51.41 
42.97 


+2.8 
-2.5 
+2.9 
-2.4 


±4.2 
±4.2 


230 


60.87 


+1.3 
-1.1 


±4.5 


36.15 


+2.8 
-2.4 


±4.8 


240 


51.45 


+1.1 
-1.2 


±4.0 


30.46 


+2.6 
-2.4 


±4.4 


250 


43.68 


+1.1 
-1.1 


±4.0 


25.81 


+2.7 
-2.2 


±4.2 


260 
270 


37.25 
31.91 


+1.3 
-1.1 
+1.2 
-1.3 


±4.0 
±3.8 


21.94 
18.70 


+2.6 
-2.3 
+2.4 
-2.2 


±4.5 
±4.3 


280 
290 
300 


27.39 
23.65 
20.47 


+0.9 
-1.4 
+1.2 
-1.2 
+1.1 
-1.3 


±4.4 
±4.2 
±4.5 


16.02 
13.79 
11.90 


+2.4 
-1.9 
+2.3 
-2.0 
+2.2 
-1.8 


±4.9 
±4.5 
±5.0 



Table 1: Numerical values for the total 
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cross section at 



LHC at 7TeV. 



Mh 

[GcV] 



a{WH) 

[fb] 



Pert 



PDF+q;^ 

[%] 



a{ZH) 



Pert 



100 
105 
110 
115 
120 
125 
130 
135 
140 
145 
150 
155 
160 
165 
170 
175 
180 
185 
190 
195 
200 
210 
220 
230 
240 
250 
260 
270 
280 
290 
300 



3035 
2625 
2273 
1976 
1731 
1523 
1342 
1183 
1048 
933.1 
827.5 
736.3 
644.0 
594.0 
534.3 
483.9 
434.5 
402.8 
366.0 
334.6 
305.5 
257.0 
217.5 
185.9 
158.9 
136.8 
118.3 
102.8 
89.49 
78.61 
68.86 



+1.3 

-0.9 
+1.1 
-0.8 
+0.8 
-1.1 
+1.2 
-0.6 
+1.1 
-0.7 
+0.9 
-1.0 
+1.0 
-0.8 
+1.2 
-0.9 
+0.8 
-1.1 
+1.1 
-0.8 
+0.9 
-1.1 
+0.9 
-1.0 
+0.8 
-0.9 
+0.8 
-1.0 
+0.9 
-1.1 
+ 1.1 
-0.8 
+1.1 
-1.0 
+1.0 

-1.1 

+0.9 
-1.1 
+0.9 
-1.2 

+1.0 
-1.0 
+0.9 
-1.1 
+1.3 
-1.0 
+1.1 
-0.9 
+1.1 
-1.0 
+0.8 
-1.2 
+0.9 
-1.2 
+1.2 
-1.1 
+ 1.1 
-1.1 
+ 1.0 
-1.1 
+1.3 
-1.1 



±3.7 
±3.5 
±3.8 
±3.8 
±3.8 
±3.8 
±3.3 
±2.9 
±3.1 
±3.3 
±2.7 
±3.1 
±3.1 
±2.4 
±2.8 
±2.9 
±2.8 
±2.5 
±2.8 
±2.7 
±3.0 
±2.6 
±2.8 
±3.5 
±3.3 
±3.0 
±2.8 
±2.6 
±3.0 
±3.2 
±3.3 



1700 
1483 
1297 
1142 
1008 
893.2 
793.9 
706.6 
633.4 
567.2 
508.2 
457.3 
404.1 
375.6 
340.2 
308.8 
278.5 
256.2 
233.6 
214.4 
196.6 
165.4 
140.3 
119.3 
101.7 
86.96 
74.80 
64.48 
55.83 
48.67 
42.43 



+2.3 
-1.9 
+2.2 
-2.1 
+2.6 
-2.0 
+2.9 
-1.9 
+2.9 
-2.2 
+3.2 
-2.2 
+3.4 
-2.2 
+3.4 
-2.6 
+3.4 
-2.6 
+3.8 
-2.5 
+3.8 
-2.4 
+3.8 
-2.8 
+4.1 
-2.8 
+4.3 
-2.7 
+4.1 
-2.9 
+4.0 
-3.1 
+4.3 
-3.3 
+4.1 
-3.3 
+4.1 
-3.3 
+4.1 
-3.4 
+4.2 
-3.5 
+4.5 
-2.9 
+4.1 
-3.2 
+4.1 
-3.1 
+3.9 
-3.0 
+3.7 
-2.9 
+3.8 
-2.6 
+3.3 
-2.5 
+3.3 
-2.4 
+3.0 
-2.2 
+2.9 
-2.2 



Table 2: Numerical values for the total 
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cross section at 



LHC at 14 TeV. 



Mh 


a{WH) 


Pert 


PDF+a^ 


a{ZH) 


Pert 


PDF+Qis 


[GeV] 


[fb] 


[%] 


[%] 


[fb] 


[%] 


[%] 


100 


278.0 


+4.0 
-2.9 


±5.1 


161.2 


+4.1 
-3.0 


±5.5 


105 


235.7 


+4.1 

-2.8 


±5.3 


138.2 


+4.3 
-3.0 


±5.6 


110 


200.9 


+4.0 
-2.9 


±5.5 


119.0 


+4.2 
-3.1 


±5.7 


115 


172.1 


+4.0 
-3.0 


±5.5 


103.0 


+4.2 
-3.1 


±5.8 


120 


148.0 


+4.1 
-3.0 


±5.8 


89.46 


+4.3 
-3.2 


±5.9 


125 


127.5 


+4.0 
-3.1 


±6.0 


77.86 


+4.2 
-3.2 


±6.1 


130 


110.2 


+4.0 
-2.9 


±6.2 


67.98 


+4.3 
-3.1 


±6.0 


135 


95.61 


+4.1 
-3.0 


±6.1 


59.55 


+4.4 
-3.2 


±6.0 


140 


83.06 


+4.2 
-3.1 


±6.4 


52.23 


+4.5 
-3.3 


±6.2 


145 


72.36 


+4.3 
-3.0 


±6.5 


46.01 


+4.6 
-3.2 


±6.3 


150 


63.16 


+4.4 

-3.1 


±6.3 


40.59 


+4.7 
-3.3 


±6.0 


155 


55.28 


+4.4 
-3.1 


±6.9 


35.91 


+4.8 
-3.3 


±6.6 


160 


48.49 


+4.6 

-3.0 


±6.9 


31.84 


+4.9 
-3.3 


±6.5 


165 


42.54 


+4.5 
-3.2 


±7.8 


28.23 


+4.9 
-3.4 


±7.3 


170 


37.51 


+4.6 
-3.3 


±7.1 


25.15 


+5.0 
-3.5 


±6.6 


175 


33.26 


+4.7 
-3.4 


±6.9 


22.42 


+5.2 
-3.6 


±6.4 


180 


29.53 


+4.9 
-3.4 


±7.3 


20.02 


+5.3 
-3.6 


±6.6 


185 


26.29 


+4.9 
-3.4 


±7.7 


17.92 


+5.4 
-3.7 


±6.9 


190 


23.39 


+5.0 

-3.6 


±7.9 


16.03 


+5.4 
-3.8 


±7.0 


195 


20.86 


+5.1 
-3.6 


±7.7 


14.38 


+5.6 
-3.9 


±6.8 


200 


18.61 


+5.0 
-3.6 


±7.6 


12.91 


+5.5 
-3.9 


±6.5 



Table 3: Numerical values for the total cross section at the Tevatron. 
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(a) (b) 

Figure 8: K-factor at the LHC at (a) 7TeV and (b) 14TeV center-of-mass energy, 
including NNLO QCD and electro-weak corrections, with and without the newly 
evaluated top-quark induced terms. 




(a) (b) 

Figure 9: K-factor at the Tevatron for (a) WH and (b) ZH production, including 
NNLO QCD and electro-weak corrections, with and without the newly evaluated 
top-quark induced terms. 
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